Evidence for a critical velocity in a Bose— Einstein condensed gas 
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We have studied dissipation in a Bose-Einstein condensed gas by moving a blue detuned laser beam 
through the condensate at different velocities. Strong heating was observed only above a critical 
velocity. 
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Macroscopic quantum coherence and collective excita- 
tions are key features in our understanding of the phe- 
nomenon of superfluidity. The superfluid velocity is pro- 
portional to the gradient of the phase of a macroscopic 
wavefunction. Collective excitations determine a criti- 
cal velocity below which the flow is dissipationless. This 
velocity is given by Landau's criterion 



p 



(1) 



where e is the energy of an excitation with momentum 
p. Critical velocities for the breaking of Cooper pairs in 
3 He and the generation of rotons || and vortices || in 
4 Hc have been extensively studied. 

Bose-Einstein condensed gases (BEC) are novel quan- 
tum fluids 0]. Previous work has explored some aspects 
related to superfluidity such as the macroscopic phase || 
and the phonon nature of low-lying collective excitations 
[^],[|. In this Letter we report on the measurement of 
a critical velocity for the excitation of a trapped Bose- 
Einstein condensate. In analogy with the well known 
argument by Landau and the vibrating wire experiments 
in superfluid helium W, we study dissipation when an 
object is moved through the fluid. Instead of a massive 
macroscopic object we used a blue detuned laser beam 
which repels atoms from its focus to create a moving 
boundary condition. 

The experiment was conducted in a new apparatus for 
the production of Bose-Einstein condensates of sodium 
atoms. The cooling procedure is similar to previous work 
H — the new features have been described elsewhere ||. 
Briefly, laser cooled atoms were transferred into a mag- 
netic trap in the Ioffe-Pritchard configuration and further 
cooled by rf evaporative cooling for 20 seconds, resulting 
in condensates of between 3 and 12 x 10 6 atoms. After 
the condensate was formed, we reduced the radial trap- 
ping frequency to obtain condensates which were consid- 
erably wider than the laser beam used for stirring. This 
decompression was not perfectly adiabatic, and heated 
the cloud to a final condensate fraction of about 60%. 
The final trapping frequencies were v r = 65 Hz in the 
radial and v z = 18 Hz in the axial direction. The re- 
sulting condensate was cigar-shaped with Thomas-Fermi 



diameters of 45 and 150 fim in the radial and axial direc- 
tions, respectively. The final chemical potential, transi- 
tion temperature T c and peak density uq of the conden- 
sate were 110 nK, 510 nK and 1.5 x 10 14 cm -3 , respec- 
tively. 

The laser beam for stirring the condensate had a wave- 
length of 514 nm and was focused to a Gaussian 1/e 2 
beam diameter of 2w = 13/im. The repulsive optical 
dipole force expelled the atoms from the region of high- 
est laser intensity. A laser power of 400 /iW created a 
700 nK barrier resulting in a cylindrical hole ~ 13/im in 
diameter within the condensate. The laser barrier cre- 
ated a soft boundary, since the Gaussian beam waist 
was more than 10 times wider than the healing length 
£ = (87ran ) -1 / 2 = 0.3/im, a being the two-body scatter- 
ing length. 





a) 



b) 



FIG. 1. Stirring a condensate with a blue detuned laser 
beam, a) The laser beam diameter is 13/im, while the radial 
width of the condensate is 45/xm. The aspect ratio of the 
cloud is 3.3. b) In situ absorption image of a condensate 
with the scanning hole. A 10 kHz scan rate was used for 
this image to create the time-averaged outline of the laser 
trajectory through the condensate. 

The laser was focused on the center of the cloud. Us- 
ing an acousto-optic deflector, it was scanned back and 
forth along the axial dimension of the condensate (fig. 
1). We ensured a constant beam velocity by applying 
a triangular waveform to the deflector. The beam was 
scanned over distances up to 60/im, much less than the 
axial extent of the condensate of 150/im. Therefore, ax- 
ial density variations were small within the scan range. 
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The scan frequencies were varied by a factor of 3 be- 
tween 56 and 167 Hz. Scan velocities close to the speed 
of sound required scan frequencies much larger than the 
axial trapping frequency v z . 

After exposing the atoms to the scanning laser for 900 
ms, we allowed the cloud to equilibrate for 100 ms, then 
turned off the magnetic trap and recorded the time-of- 
flight distribution after 35 ms on a CCD camera using 
near-resonant absorption imaging. The condensate frac- 
tion Nq/N was determined from fits to the bimodal ve- 
locity distribution ]l0[ |. We found that the decreasing 
condensate fraction was a more robust measure of the 
heating induced by the moving laser beam than the tem- 
perature extracted from the wings of the thermal cloud. 
No heating was observed when the laser beam was kept 
stationary for the entire 900 ms. 

Figure 2 shows the effect of the moving laser on the 
condensate for three different scan rates / = 56, 83 and 
167 Hz. The heating rate is higher (larger final thermal 
fraction) for larger drive amplitudes d and higher scan 
frequencies /. When the same data are replotted as a 
function of the velocity of the laser beam 2df (fig. 3), two 
features emerge immediately. First, all the three data 
sets collapse onto one universal curve, indicating that 
the heating of the condensate depends primarily on the 
velocity of the beam, and not on either frequency or am- 
plitude independently. This suggests that the observed 
dissipation is not strongly affected by the trapping po- 
tential and discrete resonances, but rather, reflects bulk 
properties of the condensate. 
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FIG. 2. Heating a condensate with a moving laser beam. 
The final thermal fraction increased with both the scan ampli- 
tude and the scan rate of the laser beam. The total exposure 
time was the same for all data points. They represent single 
shot measurements. Solid lines are smoothing spline fits to 
guide the eye. 

The second feature is that we can distinguish two 
regimes of heating separated by the dashed line in fig. 
3. For low velocities, the dissipation rate was low and 
the condensate appeared immune to the presence of the 
scanning laser beam. For higher velocities, the heating 
increased, until at a velocity of about 6 mm/s the conden- 



sate was almost completely depleted for a 900 ms expo- 
sure time. The cross-over between these two regimes was 
quite pronounced and occurred at a velocity of about 1.6 
mm/s. This velocity should be compared with the speed 
of sound in the condensate. Since a condensate released 
from the magnetic trap expands with a velocity propor- 
tional to the speed of sound, we could determine its value 
directly from time-of-flight absorption images to be 6.2 
mm/s (at the peak density), almost a factor of 4 larger 
than the observed critical velocity. 
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FIG. 3. Evidence for a critical velocity. The same three 
data sets as in figure 2 are replotted vs. laser beam velocity 
and all appear to lie upon a universal curve. The dashed line 
separates the regimes of low and high dissipation. The right 
axis indicates the temperature T = (1 - N Q /N) 1/3 T C . The 
peak sound velocity is marked by an arrow. The data series 
for 83 and 167 Hz showed large shot-to-shot fluctuations at 
velocities below 2 mm/sec. The solid line is a smoothing 
spline fit to the 56 Hz data set to guide the eye. 

To rule out the possibility of heating through the non- 
condensed fraction, a control experiment was performed 
on clouds at two different temperatures above T c , 800 
nK and 530 nK, the latter quite close to the transition 
temperature. No heating of the cloud was observed in 
either case, for scan velocities of up to 14 mm/s; this 
reflects the small overlap between the scanning laser and 
the large thermal cloud. Since these clouds are typically 
not in the hydrodynamic regime of collisions, we used 
a single particle model for heating based on collisions 
with a moving wall to scale our measurements from above 
to below T c . We obtained a conservative upper bound 
of 15% of the observed temperature rise which could be 
attributed to the non-condensed fraction. 

What determines the critical velocity? In a quasi- 
homogeneous dilute Bose gas near zero temperature, the 
critical velocity for phonon excitation is the speed of 
Bogoliubov sound, which depends on the density n(r) 
through 



Unh 2 a 



,(r) 



(2) 



where M is the atomic mass. Alternately, one may ex- 
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pect that the rapid flow around the laser beam generates 
vorticity in the fluid. For the excitation of vortex pairs 
in a channel of diameter D, eq. (1) leads to a critical 
velocity Q 



h (D 



(3) 



Since this velocity scales inversely with the size of the 
system, it is lower than the value predicted by eq. 2. For 
our situation eq. 3 can only give an approximate esti- 
mate due to the inhomogeneous density and the different 
geometry. For a channel width of 15/im and a healing 
length of 0.3/^m it yields a critical velocity of about 0.7 
mm/s, a factor of 2 below our measurement. 

The criterion given in eq. 3 is only based on consid- 
erations of energy and momentum and provides a lower 
bound to v c . In addition, one must consider how to pro- 
duce the excitations dynamically. Several mechanisms 
for creating vortex lines have been discussed, including 
remanent vorticity and vorticity pinned to the surface 
[ p2[ . In the latter case, the critical velocity depends on 
the surface roughness. The gaseous condensates are con- 
fined to magnetic traps which provide perfectly smooth 
boundary conditions. Therefore, we expect that the crit- 
ical velocity is determined by nuclcation dynamics rather 
than by purely energetic arguments. 

The relevant criterion for the onset of dissipation in 
a Bose condensed gas obeying the nonlinear Schrodinger 
equation has been discussed in several papers. Accord- 
ing to these theories, dissipation ensues when the rel- 
ative velocity between the object and the fluid exceeds 
the speed of sound cb locally [p3| P~5|| - For an incompress- 
ible flow around a cylindrical object this velocity peaks at 
the side, reaching twice the object's speed. The hydrody- 
namic equations for the compressible condensate imply 
that the faster the velocity of the flow field, the lower 
is the condensate density. This effect lowers the critical 
velocity for a cylindrical object even further fll3||: 



= 0.42c l; 



(4) 



This result is independent of the size of the object and 
was corroborated by numerical simulations of the nonlin- 
ear Schrodinger equation in a homogeneous gas [p^|-p^[. 
For our conditions, this estimate yields a critical velocity 
of 2.6 mm/s, or 1.6 times the observed threshold. How- 
ever, the finite size of the condensate in our experiments 
p6| , its inhomogeneous density distribution and the soft 
boundary of the laser beam are not accounted for by the 
theory. All these effects should lower the critical velocity. 

What happens to the condensate above v c ? Numer- 
ical simulations of the nonlinear Schrodinger equation 
were used to study the flow field around an object mov- 
ing through a homogeneous condensate |I^,|l4|,[l7j . These 
studies show that above a critical velocity given by Eq. 



(4) the superfluid flow becomes unstable against the for- 
mation of quantized vortex lines, which signals the onset 
of a new, dissipative regime. Pairs of vortices with op- 
posite circulation are generated at opposite sides of the 
object in order to reduce the high local flow speed. The 
rate of heating can be estimated from the energy of vor- 
tices and the vortex shedding frequency. 

The energy of a vortex pair e vp is estimated with the 
assumption that the vortices are separated by the object 
diameter 2w [fill, 



„nh 2 2w 
e vp = 2nD^-j- In — 



(5) 



where D is the radial width of the condensate and n is 
the density. This yields a value of about 3.4 mK or 1.3 
nK/atom. Numerical simulations have shown that 

for v > v c the rate of vortex pair shedding is proportional 
to!)-D c . The proportionality constant, besides a numer- 
ical factor, is the mean field energy in frequency units di- 
vided by the speed of sound. Thus the rate of change in 
temperature should have the form T = k(v — v c ). Using 
the estimate for the vortex energy, the model predicts 
k w 160 nK/mm, in rough agreement with our mea- 
sured heating rate near the threshold which gives k = 62 
nK/mm. 

In conclusion, we have established a method for study- 
ing dissipation in a Bose condensate by implementing a 
scanning "hole" induced by a far-off-resonant laser beam. 
Both the laser beam and the trapped condensate provide 
clean boundary conditions. This and the simplicity of 
the system makes it amenable to theoretical treatments 
of vortex nucleation and dissipative dynamics. We found 
evidence for a critical velocity for excitation of the con- 
densate. Both the onset and the magnitude of the ob- 
served heating are in qualitative agreement with model 
calculations based on the non-linear Schrodinger equa- 
tion which predict dissipation when the flow field be- 
comes locally supersonic. In contrast, a similar study 
on rotational stirring concluded that the speed of sound 
near the stirrer is irrelevant for vortex nucleation fl9| . In 
further studies, we plan to vary the geometry and density 
in order to distinguish between different predictions for 
critical velocities which depend only on geometry (eq. 3), 
or only on the density (eq. 4). All the calculations were 
done in two dimensions and at zero temperature. In the 
experiment, the laser beam passes through the surface 
of the condensate where the density vanishes. Because 
of this and the non-zero temperature, we expect finite 
dissipation even at low velocities and a smooth crossover 
between low and high dissipation. More precise mea- 
surements of the heating should allow us to study these 
finite-size and finite-temperature effects. 
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